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November  1955 


EFFECT  OF  A  HEMISPHERICAL  BASE  ON  THE  AERODYNAMIC  CHARACTERISTICS  OF  SHELL 

ABSTRACT 

The  addition  of  a  hemispherical  base  to  a  square  based  model 
produces  marked  dynamic  instability.  A  comparison  of  the  aerodynamic 
characteristics  of  square  and  hemispherical  based  models  is  given. 

Estimates  of  the  damping  and  Magnus  forces  of  the  hemispherical  based 
model  are  obtained.  The  reason  for  the  instability  is  quite  clear 
although  the  mechanism  uhich  produces  this  result  seems  quite  complex. 
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TABIE  OP  SYMBOLS 

Center  of  mass  of  projectile 

Center  of  pressure  of  normal  force 

Center  of  pressure  of  the  increment  of  force  between 
the  square  and  hemispherical  based  models 

Diameter  of  projectile 

Drag  force  coefficient 

IC  when  6  =  0 


Kdb  /  6  =  0 

Lift  force  coefficient 

Normal  force  coefficient 

Overturning  moment  coefficient  f 

Magnus  force  coefficient 

Magnus  moment  coefficient 

Damping  force  coefficient 

Damping  moment  coefficient 

Spin  deceleration  moment  coefficient 

Square  based  model  coefficient  corrected  to  a  location 
corresponding  to  the  c.m.  of  the  equivalent  length 
hemispherical  based  model. 

Estimated  coefficient  of  the  hemispherical  baaed  model 

Axial  radius  of  gyration 

Transverse  radius  of  gyration 

Amplitude  of  fast  rate  of  yawing  motion 

Amplitude  of  slow  rate  of  yawing  motion 

Mass  of  projectile 

Mach  number 


.2  2 

A 


4fiPd  «i 


gyroscopic  stability  factor 


i 


C to 


62 


2(K  -  k^2  K_) 

- 5 - 7j —  dynamic  stability  factor 

*L+V  KA 

Radius  of  lift  swerving  motion 

Exponential  damping  coefficient  of  fast  rate 
Exponential  damping  coefficient  of  slow  rate 
Magnitude  of  yaw  single 
Mean  yaw  squared 

Standard  error  for  fit  of  yaw  curve  to  data 
Standard  error  for  fit  of  swerving  motion  curve  to  data 
Standard  sea  level  air  density 
Mr  density 
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INTRODUCTION 


* 

A  program  of  20nm  models  with  square  and  hemispherical  bases  was  fired 
in  Exterior  Ballistics  Laboratory' s  small  aerodynamic  free  flight  spark 
range  in  order  to  try  to  determine  comparative  aerodynamics  of  these 
configurations .  The  results  obtained  from  the  32  rounds  which  could  be 
reduced  are  given  in  this  report.  The  four  types  of  models  fired  con¬ 
sisted  of  the  basic  5 <183  calibers  long,  square  based  model  and  three 
modifications:  (l)  the  addition  of  a  hemispherical  base,  (2)  the 
addition  of  a  0.556  caliber  cylinder,  and  (3)  the  addition  of  a  0.556 
caliber  cylinder  plus  a  hemispherical  base.  The  second  modification 
was  made  to  determine  whether  the  added  length  of  the  hemisphere  was 
the  cause  of  instability,  and  the  third  modification  was  made  to  sea 
if  the  addition  of  the  hemisphere  to  the  longer  model  gave  the  same 
type  of  changes  as  the  first  modification.  All  four  models  can  be 
seen  in  Figure  1  with  the  dimensions  being  given  in  Figure  2.  The 
customary  methods  of  data  reduction^  were  used. 

AERODYNAMIC  COEFFICIENTS 


The  drag  coefficient  as  obtained  from  the  least  squares  fit  of  a 
cubic  equation  in  distance  down  range  to  the  time  interval  contains  the 
effect  of  the  variation  of  this  coefficient  with  yaw  as  shown  by 

KD-KDo*\a'52- 

Since  there  were  very  few  rounds  which  had  different  yaws  at  essentially 
the  same  Mach  number,  it  was  not  possible  to  determine  2  by  a  fit  of 

-  2  5 

Kp  to  5  .  A  value  of  2  =  2.0  per  radian  squared  was  therefore 

6 

assumed  in  order  to  determine  the  values  of  (see  Figure  7).  Values 
of  the  spin  deceleration  coefficient,  K,,  could  not  be  determined 
because  pins  for  a  spin  reduction  were  not  placed  in  the  bases  of  the 
models . 

v. 

Additional  data  on  hemispherical  based  shell  can  be  found  in  Reference  7- 
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Although  there  is  a  variation  of  with  the  yaw,  the  yaw  for  this  set 
of  rounds  is  so  small  that  the  correction  has  been  neglected.  The  curves 
for  Kjj  as  shown  in  Figure  8  seem  to  be  quite  well  defined  and  show  an 
appreciable  difference  between  the  square  and  hemispherical  based  models 
which  is  not  attributable  to  experimental  inaccuracy.  It  was  felt  that 
the  values  of  for  corresponding  bases  were  within  the  accuracy  of 
determination,  therefore  only  one  curve  for  each  type  of  base  was  drawn. 
There  is  not  enough  data  in  the  transonic  region  to  accurately  deter¬ 
mine  the  shape  of  the  or  the  curve  in  the  neighborhood  of  M  -  1. 

The  general  trend  of  the  curves  of  and  Kj^  vs.  Mach  number  have, 
therefore,  been  drawn  to  agree  with  Reference  2.  Since  the  overturning 
moment  coefficient  is  contaminated  by  the  different  c.m.  locations  for 
the  different  models,  the  values  of  for  the  square  based  models  have 
been  evaluated  at  a  point  which  is  the  same  distance  from  the  nose  of 
the  c.m.  of  the  hemispherical,  based  model  of  comparable  length  and  are 
denoted  (see  Figure  9) •  The  plot  of  the  center  of  pressure  of  the 
normal  force  shown  in  Figure  10,  together  with  K^,  are  probably  even 
more  descriptive  of  the  effect  of  the  hemisph  arlcal  base  on  the  over¬ 
turning  moment  than  the  plot  of  K^. 


The  curves  of  the  Magnus  moment  in  Figure  11  are  some  of  the  better 
determined  curves  obtained  from  this  firing,  Again  there  is  a  distinct 
difference  in  the  values  of  the  hemispherical  and  square  based  models 
especially  in  the  range  of  M  y  1.1.  Again  it  was  felt  that  only  one 
curve  for  each  base  type  was  warranted.  Since  the  swerve  due  to  the 
Magnus  force  for  all  of  these  rounds,  was  less  than  the  accuracy  of 
measurement  the  determinatiu i  of  Kp  with  any  degree  of  accuracy  was 
impossible . 


The  damping  moment  coefficients  in  Figure  12  also  show  only  a 
trend,  but  they  do  indicate  a  distinct  difference  in  the  two  types  of 
bases  which  is  not  attributable  to  experimental  error.  Values  of  the 
damping  force  could  not  be  obtained  since  models  of  the  same  configura¬ 
tion  but  with  different  center  of  mass  location  were  not  fired. 


8 


The  values  of  the  coefficients  near  Mach  =  1  may  be  affected  by  the 
interference  of  reflected  nose  shock  waves  with  the  afterbody  of  the 
model,  therefore  there  is  some  doubt  as  to  their  accuracy. 

DISCUSSION 

The  outstanding  difference  in  the  shadowgraphs  of  the  square  and 
hemispherical  based  models  is  the  flow  over  part  of  the  hemisphere  and 
the  resulting  shock  wave  when  the  boundary  layer  separates  from  the 
base,  as  is  shown  in  Figures  5  and  4. 

There  is  a  very  distinct  difference  in  the  drag  coefficients  of 
the  square  and  hemispherical  based  models,  with  the  hemispherical  base 
displaying  a  noticeably  greater  drag  than  the  square  based  models 
(see  Figure  7).  The  drag  for  a  short,  large  angled  boattail  is  higher 
than  that  for  a  square  based  missile  as  is  shown  in  Reference.  4  and  is 
theoretically  discussed  by  J.  Sternberg^ \  The  separation  of  the  flow 
from  the  hemispherical  base  in  essence  gives  a  short,  large  angled 
boattail. 

The  normal  force  is  increased  by  the  addition  of  a  hemispherical 
base  to  the  rear  of  a  regular  square  base.  The  largest  difference  in 
Kjj  for  the  two  types  of  bases  is  in  the  transonic  region.  Since  the 
center  of  pressure  of  the  normal  force  is  moved  rearward  by  the 
addition  of  the  hemispherical  base  as  seen  in  Figure  10,  it  seems  that 
there  has  been  an  increase  in  the  pressure  difference  at  the  rear  of 
the  projectile.  This  pressure  difference  could  be  the  result  of  the 
difference  in  separation  points  of  the  boundary  layer  on  opposite 
sides  of  the  hemispherical  base  and  the  resulting  shock  wave.  The 
separation  angles,  which  are  defined  as  the  angles  made  by  the  radius 
vector  from  the  center  of  the  sphere  to  the  point  of  separation  and 
the  model’s  axis  as  illustrated  in  Figures  5  and  6,  were  measured  on 
numerous  plates .  Only  in  a  very  few  cases  measurable  differences  in 
the  separation  angles  on  the  two  sides  of  the  projectile  were  obtained. 
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Two  photographs  illustrating  a  large  difference  in  separation  angles 
are  shown  in  Figures  5  and  6  where  a  difference  of  3°  and  2°  respectively, 
was  measured. 


The  largest  difference  in  the  normal  forces  and  their  centers  of 
pressure  as  well  as  in  the  overturning  moments  is  in  the  transonic  region 
with  the  difference  decreasing  for  increasing  Mach  number.  This  indicates 
that  there  are  changes  in  the  flow  forward  of  the  base  due  to  the 
additional  shock  wave  at  the  point  of  the  boundary  layer  separation  from 
the  hemispherical  base.  This  disturbance  could  be  propagated  forward 
in  the  boundary  layer  and  also  in  the  potential  flow  at  transonic  Mach 
numbers.  With  the  aid  of  this  observation  the  change  in  the  normal 
force  coefficient  due  to  the  addition  of  a  hemispherical  base  can  be 
checked  by  means  of  the  change  in  the  overturning  moment  coefficient 
and  some  assumptions  about  the  location  of  the  normal  force. 


Since  the  overturning  moment  coefficients  for  the  square  based 
models  have  been  evaluated  at  the  position  of  the  c.m.  of  the  corresponding 
hemispherical  based  models,  the  difference  in  the  coefficients  should 
be  completely  due  to  the  change  in  the  flow  caused  by  the  addition  of 
the  boat  tail.  Monk's  linearized  slender  body  theory  predicts  a  change 
in  the  normal  force  on  a  projectile  corresponding  to  the  change  in  its 
cross  sectional  area.  Therefore,  it  seems  reasonable  to  assume  that 
the  coefficient,  which  represents  the  change  in  normal  force 

coefficient  due  to  the  addition  of  the  hemisphere  will  act  at  a  point 
in  the  region  in  which  the  projectile  is  undergoing  the  decrease  in 
cross  sectional  area.  The  effective  boattall  is  from  the  beginning  of 
the  hemisphere  to  the  separation  point,  and  is  about  0.10  calibers  in 
length.  It  will  be  assumed  that  acts  at  the  middle  of  this  region. 

For  the  short  hemispherical  based  model,  where  the  beginning  of  the 
hemisphere  is  1.63  calibers  behind  the  c.m.,  the  distance  from  the 
c.m.  to  £K_  is 


=  1.70  +  0.2  calibers. 
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iWiv-pi-r. 


and  for  the  long  hemispherical  based  model,  where  the  beginning  of  the 
hemisphere  is  1.92  calibers  behind  the  c.m. 


1.97  +  0-2  calibers. 


where  the  +  0.2  caliber  is  used  to  give  a  probable  limit  to  the  value. 


If  the  value  of  the  change  in  between  the  square  and  hemispherical 
based  models  of  corresponding  lengths  is  obtained  from  the  graph  in  the 
appendix  and  it  is  assumed  that 


where  K.  is  the  moment  coefficient  for  the  hemispherical  based  model, 

®  r* 

then  AKg  can  be  computed.  Using  this  AK^  and  +  AKg,  the  following 

table  shows  a  comparison  of  predicted  and  observed  values  of  ICg. 


Model 

Mach 

Calculated 

Calculated 

Observed 

Typ«* 

Humber 

*H 

** 

ss 

0.9 

0.18  +  0.02 

1.11  + 

0.02 

1.12 

IS 

0.9 

0.18  +  0.02 

1.11  + 

0.02 

1.12 

SH 

1.6 

-0.006  +  0.001 

1.014  + 

0.001 

1.10 

IS 

1.6 

0.000 

1.02 

1.10 

These  show  good  agmement  with  the  observed  values ,  especially  at 


M  -  0.9. 


The  damping  moment  coefficient  for  the  hemispherical  based  models  is 
markedly  negative  in  the  transonic  region  and  seems  to  be  approaching  a 
zero  or  positive  value  as  the  Mach  number  increases  as  shown  in  Figure  12. 
Since  a  positive  Kg  indicates  that  there  is  a  resistance  to  the  change 
in  yaw,  the  negative  Kg  means  'hat  the  amount  of  yaw  is  being  increased. 
Under  these  circumstances  the  projectile  is  unstable. 

On  the  basis  of  the  good  agreement  obtained  for  the  estimate  of 
the  normal  force  coefficient,  it  is  felt  that  a  fairly  good  estimate  of 
Kg  can  also  be  obtained  in  the  same  manner. 

*  Model  types  are  identified  in  Table  I  of  the  Appendix. 
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The  damping  force  is  considered  to  be  the  6ame  type  of  phenomenon 
as  the  normal  force,  i.e.,  caused  by  the  change  in  the  momentum  of  the 
potential  flow.  The  velocity  field,  however,  is  the  result  of  the 
velocity  induced  by  the  cross  spin  rather  than  the  cross  velocity  due 
to  yaw.  It  therefore  seems  reasonable  to  assume  that  the  location  of 
AKg  is  at  the  same  place  as  Although  Kg  for  the  square  based 

models  is  not  known,  a  comparison  of  the  known  force  coefficients  of 
this  report  with  those  in  References  2  and  3  indicate  that  the  forces 
for  the  two  slightly  different  models  are  about  the  same  within  about 
5*. 


Therefore,  using  References  2  and  3,  it  is  assumed  that  Kg  =  -  6.0 
for  the  short  and  long  square  based  models  at  M  *  0.9  •  At  M  »  1.6  this 

based  model  K^  values 
of  this  report,  which  are  at  the  c.m.  of  the  square  based  models,  cure 
evaluated  at  a  point  corresponding  to  the  c.m.  of  the  hemispherical 
based  models  and  listed  in  the  following  table. 


value  is  -  4.0.  These  Kg  values  and  the  square 


Model 

Mach 

V 

V 

Type 

Number 

SS  and  U5 

0.9 

-  5.8 

-  2.6 

SS  and  I£ 

1.6 

-  3.8 

4.2 

Similar  to  the  previous  discussion  it  is  assumed  that: 

V  V  -  (c'p-acs)a's 

and 

=  V  + 

The  following  table  gives  computed  values  of  1L. 


Madel 

Mach 

Calculated 

Calculated 

Type 

Number 

S 

SH  and  LH 

0.9 

2.5  +  0.5 

-  5.5  +  0.! 

SH  and 

1.6 

3.1  +  0.6 

-  0.7  +  0.< 

These  estimates  Indicate  that  the  change  in  damping  force  in  a  positive 
direction  is  quite  pronounced. 


Figure  11  shows  that  the  addition  of  the  hemispherical  base  causes 
the  Magnus  moment  to  assume  a  large  positive  value  even  in  the  supersonic 
region.  This  condition  alone  is  sufficient  for  instability. 

In  Reference  6  it  can  be  seen  that  the  limits  for  dynamic  stability 

are 

(1)  ^  +  k2-2  Kh  -  if2  Ka7  o, 

(2)  04  6  4  2, 

(5)  6  *  i(2  -V  ' 

where  s  is  the  gyroscopic  stability  factor  and  s'  is  the  dynamic  stability 
factor.* 


As  shown  in  Table  III  of  the  Appendix,  all  of  the  hemispherical 
based  models  are  gyros copically  stable  (s  7  l).  Only  two  of  the  models, 
however,  satisfy  the  dynamic  stability  condition  that  0  <  ¥  <  2,  and 
these  two  models  do  not  satisfy  the  condition  that  s -  b(2  -  s) — * 

It  is  interesting  to  note  that  even  if  were  positive  for  the  hemis¬ 
pherical  based  models,  they  would  still  not  be  stable  because  of  the 
large  positive  Magnus  moment.  Average  values  of  the  coefficients  for 
the  hemispherical  based  models  for  this  report  are  as  follows: 


*L*l-° 

~0*5 

KjjSjO.15 

k^2  ^  10 

K^O.4 

kg"2~:0.6 

Kg^-4.0 

ka  a*  0.01 

With  the  aid  of  the  conditions  (l)  and  (2)  for  dynamic  stability, 
an  interesting  observation  can  be  made.  If  Kg  =  4.0  instead  of  being 
negative  then 

1^  +  k2‘2  -  kx"2  Ka  1.0  +  0.6  (+4.0)  -  10(0.01)  «  5-5 


♦Algebraic  definition  of  s  is  given  in  the  Table  of  Symbols. 
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and  the  first  condition  is  satisfied.  The  dynamic  stability  factor 


s 


2(kl  -  y%> 

ki'"  ka 


2 


-  10(0.4)1 
3.3 


=  -  1.8 


and  the  missile  would  still  he  unstable,  since  s  is  not  between  0  and  2. 


A  further  examination  of  the  Magnus  effects  shows  an  interesting 
feature  about  them.  The  Magnus  effects  are  usually  thought  to  be  a 
boundary  layer  phenomenon.  Since  the  effects  of  the  additional  shock 
wave  on  the  hemispherical  base  would  bs  felt  upstream  in  the  subsonic 
boundary  layer  and  result  in  a  change  in  the  boundary  layer  character¬ 
istics,  it  is  conceivable  that  there  would  be  changes  in  the  Magnus 
force  and  moment.  The  change  in  the  Magnus  moment  is  verified  by 
Figure  11. 


The  Magnus  force  coefficient  is  also  estimated  by  the  same  procedure 
as  was  used  f rr  and 

ZsKg,  is  assumed  to  be  different  than  that  for  the  other  two  forces.  The 
relatively  large  rotating  band  is  oelieved  to  be  a  natural  boundary; 
therefore,  the  effect  of  the  additional  shock  wave  is  assumed  to  influence 
the  boundary  layer  flow  only  behind  the  rotating  band.  If  is 
assumed  to  act  at  the  center  of  the  region  between  the  rear  of  the 
rotating  band  and  the  point  of  boundary  layer  separation  then, 

C.?.Alr  =  1.40  +  0.2  calibers 

for  the  short  hemispherical  based  model  and 

2.P. A„  =  1.40  +  0.4  calibers 

for  the  long  hemispherical  based  model.  The  tolerance  allows  the  force 
to  have  limits  covering  almost  the  whole  effected  region.  Again  using 
References  2  and  5  it  is  assumed  that  Kp  =  0.08  and  0.15  for  M  =  0.9 
and  1.6  respectively  for  both  the  short  and  long  models  since  there  is 

/  7\ 

no  effect  of  length  for  the  length  of  models  used'1  .  The  graph  values 
of  for  the  square  based  models  are  evaluated  at  a  point  corresponding 
to  the  c.m.  of  the  hemispherical  based  models  to  give  the  following  table 


Kgj  however,  the  point  of  application  of  the  force 
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Model 

Mach 

“r 

V 

Type 

Number 

SS  and  L3 

0.9 

0.39 

0.40 

SS  and  I B 

1.6 

-0.13 

-0.10 

As  before  it  is  assumed  that 


and 


-  V  ■  <a-p-a£F>a!P 


It  is  then  found  that  the 


estimated  values  are  as  shovn  in  the  following 


table 


Model 

Type 


Mach 

Number 


SH  and  LH 
SH  and  LH 


0.9  -  0.09  +  0.02  -  0.01  +  0.02 

1.6  -0.36+0.06  -0.21+0.06 


The  change  in  K^,  seems  to  be  a  very  outstanding  indication  of  the  affect 
of  the  hemispherical  base  since  it  changes  from  a  positive  quantity  for 
the  square  based  models  to  a  negative  quantity  for  the  hemispherical 
based  models.  This  change  is  very  pronounced  at  M  =  1.6. 
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